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^ I The deconfinement transition in SU(2) gauge theory and the magnetization transition 
in the Ising model belong to the same universality class. The critical behaviour of the 

^ ■ Ising model can be characterized either as spontaneous breaking of the Z2 symmetry 

CN . of spin states or as percolation of appropriately defined spin clusters. We show that 

^ ! deconfinement in SU(2) gauge theory can be specified as percolation of Polyakov loop 

\ clusters with Fortuin-Kasteleyn bond weights, leading to the same critical exponents as 

^ ■ the conventional order-disorder description based on the Polykov loop expectation value. 

O ; 

; 1. INTRODUCTION 

p I The study of critical phenomena has always been one of the most challenging and fas- 
^ ■ cinating topics in physics. One can give many examples of systems which undergo phase 
^ ■ transitions, from familiar cases like the boiling of water in a kettle or the paramagnetic- 
• r-H . ferromagnetic transition of iron at the Curie temperature, to the more complicated case of 
rS \ the transition from hadronic matter to quark-gluon plasma which is likely to be obtained 
\ by high-energy heavy-ion experiments in the coming years. Particularly interesting are 
the second-order phase transitions, characterized by a continuous variation of the order 
parameter and a divergent correlation length at the threshold. Already in the 70's one 
tried to use percolation theory in order to give a geometrical description of the critical 
behaviour of dynamical systems undergoing second-order phase transitions Perco- 
lation has in fact several features in common with such systems: power law behaviour of 
the variables at criticality, scaling relations, universality. 

The critical behaviour of the Ising model can indeed be reformulated in terms of perco- 
lation theory: magnetization sets in when suitably defined clusters of parallel spins reach 
the dimensions of the system . In particular, the critical exponents for percolation then 
become equal to the Ising exponents. 

The paramagnetic-ferromagnetic transition of the Ising model is strictly related to the 
confinement-deconfinement transition in finite temperature SU{2) gauge theory. This was 
conjectured on the basis of effective theories and confirmed by lattice studies |0. The 
order parameter for SU{2) is the lattice average of the Polyakov loop, and it behaves like 
the magnetization in the Ising model; in particular, the critical exponents are the same. 
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These analogies are the basis of this work. We have tried to see whether it is possible 
to find a description of critical behaviour in terms of percolation also for the deconfine- 
ment transition in SU{2) gauge theory. We show that in a lattice regularization which 
corresponds to the strong coupling limit, both in two and in three space dimensions, the 
percolation of Polyakov loop clusters (taken to be suitably defined areas of Polyakov loops 
L of the same sign) leads to the correct deconfinement temperature and to the correct 
critical exponents for the deconfinement. We have achieved this result by means of com- 
puter simulations of finite temperature SU{2) gauge theory and using standard finite size 
scaling tecniques to extrapolate the results to the infinite volume limit. 

2. PERCOLATION THEORY AND THE ISING MODEL 

The percolation problem is easy to formulate: just place randomly pawns on a 
chessboard. Regions of adjacent pawns form clusters. Percolation theory deals with the 
properties of these clusters when the chessboard is infinitely large. If one of the clusters 
spans the chessboard from one side to the opposite one, we say that the cluster percolates. 

Quantitatively, one counts how many pawns belong to each cluster and calculates two 
quantities: 

• The average cluster size S, defined as: 

Here Ug is the number of clusters of size s and the sums exclude the percolating cluster; 
this number indicates how big on average the clusters are which do not percolate. 

• The percolation strength P, defined as: 

p size of the percolating cluster , . 

no. of lattice sites 

By varying the density of our pawns, a kind of phase transition occurs. We pass 
from a phase of non-percolation to a phase in which one of the clusters percolates. The 
percolation strength P is the order parameter of this transition: it is zero in the non- 
percolation phase and is different from zero in the percolation phase. 

If we want to study percolation in dynamical systems, for example the Ising model, first 
of all we must define the rule to build up the clusters (for instance we can join together all 
nearest-neighbour spins of the same sign). Then one has to find out at which temperature 
we have a spanning cluster. Let us call this temperature Tp. It turns out that for T^iTp 
the percolation variables P and S behave in the following way: 

P^{Tp-Tf^, ^~|T-Tp|-> (3) 

A. Coniglio and W. Klein demonstrated that for some special definition of cluster, 
the onset of percolation coincides with the one of magnetization; besides, Pp and 7p 
coincide respectively with the magnetization exponent /3 and the susceptibility exponent 
7. Such cluster definition had already been used by Fortuin and Kasteleyn to show that the 
partition function of the Ising model can be rewritten in purely geometrical terms as a sum 
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over clusters configurations [|10[- According to the Fortuin-Kasteleyn prescription, two 
nearest-neighbouring spins of the same sign belong to the same cluster with a probability 
p = 1 — exp{—2j3) {(3 = J/kT, J is the Ising coupling, T the temperature). The result of 
Coniglio and Klein is valid in any space dimension d {d>2) and it is independent of the 
lattice geometry (square, triangular, honeycomb, etc.) as long as it is homogeneous. 

3. POLYAKOV LOOP PERCOLATION IN SU(2) GAUGE THEORY 

Finite temperature SU{N) gauge theories describe the interactions of systems of mu- 
tually interacting gluons in thermal equilibrium (A^ is the number of colours). Systems 
containing quarks together with gluons are certainly more appealing; such a scenario will 
be explored experimentally up to very high temperatures by means of high-energy heavy 
ion collisions. Nevertheless, SU{N) gauge theories are of theoretical interest; already at 
this level one has a transition between confinement to deconfinement, going from a phase 
of bound state of gluons (glueballs) to a phase of free gluons. The order parameter of 
the confinement-deconfinement phase transition of SU{N) gauge theories is the lattice 
average of the Polyakov loop < L >. It is related to the effective potential V{r) of a 
static (mass oo) quark- ant iquark pair put in the gluons' medium at temperature T at 
a distance r from each other, when r is very big: 

I < L > I ~ lim e r ^ (4) 
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The Polyakov loop is zero in the confined phase {y{r) —>■ oo) and different from zero in 
the deconfined phase {V{r) finite). Some time ago it was conjectured by B. Svetitsky and 
L. G. Yaffe that the critical behaviour of SU{N) gauge theories has a strong relationship 
with the one of simple Z{N) spin models, with which they share a common global symme- 
try . In particular, in case of second order phase transitions, both models would belong 
to the same universality class, that is they would have the same set of critical exponents. 
One simple test of the Svetitsky- Yaffe conjecture is provided by the SU{2) gauge theory: 
numerical simulations showed that its critical exponents indeed coincide with the ones of 
the Ising model 0. 

In some respect the lattice configurations of SU{2) are similar to the ones of the Ising 
model. Instead of having a spin variable at each lattice site we have the value of the 
Polyakov loop at that site. Our aim is to see whether it is possible to build clusters of 
Polyakov loops such that their percolation indices (threshold, exponents) coincide with 
the thermal ones. In order to do that, we have to face two difficulties: 

i) the Polyakov loop is not a two-valued variable like the spin in the Ising model but 
a continuous one; its values vary in a range that, with the normalization convention we 
use, is [—1, 1]; 

ii) the SU{2) Lagrangian is not directly a function of the Polyakov loop, therefore it 
is not possible to extract from it the expression of the bond probability that we need to 
build the clusters like in the Ising model. 



In a recent work ||TT| it was proved that the point i) is not a problem. If we take an Ising 
model with continuous instead of two-valued spins, the definition of Fortuin-Kasteleyn 
clusters can be generalized by taking as a bond probability between two positive (or 
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negative) nearest-neighbour spins Si and sj the following expression: 

p{i,j) = 1 - exp{-2(3siSj) (5) 

The problem ii), however, is hard to overcome. In particular, it seems that there is 
no way to express the lattice Lagrangian of SU{2) in terms of Polyakov loops for all 
lattice regularizations. Therefore we were forced to investigate SU{2) in a special lattice 
regularization, which corresponds to the so called strong coupling limit. In this case it 
was shown that the partition function of SU{2) can be written in a form which, 
apart from a factor which depends on the group measure, is the partition function of the 
continuous Ising model studied in [|lll] . The coupling k of the spin model and the coupling 



(3 of SU (2) are related to each other by this relation: 

K ~ (/3/2)^' (6) 

{Nt fixes the temporal lattice regularization). If SU{2) is approximately a continuous 
Ising model, we can try to use the general definition of clusters of |TT]]. Therefore we have 



defined our clusters as regions of like-sign Polyakov loops connected by bonds distributed 
according to the bond weight 

Pihj) = 1 - exp{-2KLiLj), (7) 

{Li and Lj are the Polyakov loop values at the sites i and j). For Nf = 2 the approximation 
of [|12| is good and we chose to investigate numerically this special case in two and three 
space dimensions. 

4. RESULTS 

We began our percolation studies performing some test runs for different lattice sizes 
to check the behaviour of our percolation variables around criticality. Figure 1 shows the 
average cluster size S" for 2 -|- 1 SU{2) in correspondence of different lattice sizes. To get 
the critical point of the percolation transition we used the method suggested in |T^. For 



a given lattice size and a value of (3 we counted how many times we found a percolating 
cluster. This number is successively divided by the total number of configurations at 
that j3. We call this quantity percolation probability. This variable is directly a scaling 
function, analogue to the Binder cumulant in continuous thermal phase transitions. Figure 
2 shows the percolation probability for the 3-dimensional case as a function of 13 for 24^ x 2, 
30^ X 2 and 40^ x 2. 

The lines cross at the same point within the errors and that restricts further our (3 range 
for the critical threshold. Besides, since the percolation probability is a scaling function, 
we could already get clear indications about the class of critical exponents of our clusters. 
In fact, if one knows the critical point and the exponent z/, a rescaling of the percolation 
probability function of {(3 - (3cr)L^''' should give us the same function for each lattice 
size. Fig. 3 and 4 show the rescaled percolation probability for 3-dimensional SU{2) using 
j3cr = 1.8747 and two different values of the exponent z/, respectively the Ising value and 
the random percolation one. The figures show clearly a remarkable scaling for the Ising 
exponent and no scaling for the random percolation exponent. 
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Figure 1. Average cluster size as a function Figure 2. Percolation probability for 3 + 

of /3 for 2 + 1 SU{2) and different lattice 1 SU{2) as a function of f3 for three lattice 

sizes. The dashed line indicates the position sizes, 
of the physical threshold. 
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Figure 3. Rescaled percolation probability Figure 4. Rescaled percolation probability 
for 3 + 1 SU{2) using per = 1-8747 and the for 3 + 1 SU{2) using p^r = 1-8747 and the 
3- dimensional Ising exponent u = 0.629. 3-dimensional random percolation exponent 

u = 0.88. 



To evaluate the exponents f3 and 7 we performed high-statistics simulations in the range 
where the percolation probability curves cross each other. The number of measurements 
we took for each value of the coupling varies from 50000 to 100000. We used the 
method 0] to determine the values of the exponents. The final results are reported in 
Table 1 and 2. The agreement both in two and in three dimensions is good. 

5. CONCLUSIONS 

We have shown that the confinement-deconfinement phase transition in finite temper- 
ature SU{2) pure gauge theory can be described as percolation of some suitably defined 
clusters of Polyakov loops of the same sign. Our result is valid only in the strong coupling 
limit, and in order to make it general a different approach seems to be necessary. The use 
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Table 1 



Comparison of the thermal and the percolation results for 2+1 SU(2). 





Critical point 




7/z/ 


V 


L Percolation 


3.443(1) 


0.128(5) 


1.752(8) 


0.98(5) 


Spont. Symm. Breaking 


3.464(14) 


1/8 


7/4 


1 



Table 2 



Comparison of the thermal and the percolation results for 3+1 SU(2). 





Critical point 


fijv 




V 


L Percolation 


1.8747(2) 


0.528(15) 


1.985(13) 


0.632(11) 


Spont. Symm. Breaking 


1.8735(4) 


0.523(12) 


1.953(18) 


0.630(14) 


Ising Model [0] 




0.518(7) 


1.970(11) 


0.6289(8) 



of effective theories for SU (2) may help to solve the problem [[15 
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